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Thorsten Kohler 

Clarendon Laboratory, Department of Physics, University of Oxford, Oxford, 0X1 3PU, United Kingdom 

(February 1, 2008) 



O 
O 

(N 

o 
O 



> 

in 
o 



o 



o 



The stability of a Bose condensed gas with a negative 
scattering length is studied with regard to the role played by 
the initial state. A trapped ideal gas ground state is shown to 
be unstable when the Gross Pitaevskii equation still predicts 
the existence of a stable state. A possible relation to a recent 
experimental study of the critical conditions for stability or 
collapse of a Bose Einstein condensate [J. L. Roberts et al., 
Phys. Rev. Lett. 86, 4211 (2001)] is discussed. 
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The possibility of manipulating interatomic forces us- 
ing magnetic field tunable Feshbach resonances has 
provided new opportunities for the study of collision phe- 
nomena in ultra cold gases. Recent experiments with 
^^Rb Q have shown how one can tune the two body 
interaction, in a condensate from strongly repulsive to 
attractive over a wide range of different s wave scat- 
tering lengths. This makes possible the study of non 
equilibrium phenomena under a variety of controlled ini- 
tial conditions. One of the most striking experiments of 
this type is that in which a condensate, first formed with 
repulsive forces, is switched to the attractive case, and 
collapses. In an idealized sense collapse means an un- 
restricted contraction of a gas within a finite amount of 
time. As shown by Pitaevskii the time dependent 
Gross Pitaevskii equation (GPE) predicts this remark- 
able phenomenon for a trapped Bose condensed gas with 
a negative scattering length in a harmonic trap as soon as 
certain conditions, on the interaction strength, and the 
initial condensate wave function, are fulfilled. 

Very recently, the control of the interatomic force in 
^^Rb [|j has allowed an experimental determination of 
tlie critical conditions for the stability, with respect to 
collapse, of a condensate. In the course of these experi- 
ments, the interatomic force of the gas was tuned towards 
specific negative s wave scattering lengths by means of 
magnetic fields, using a Feshbach resonance. Since the 
initial scattering length was chosen to be negligibly small, 
from the point of view of the GPE the prepared initial 
state can be interpreted as an ideal gas condensate wave 
function. 

Since their first observation in experiments with ^Li , 
Bose condensed gases with a negative scattering length 
have been studied by numerous authors. Most of this 
work has used the GPE and focused on the stability con- 
ditions for a ground state solution (see, e^., |^ and ref- 
erences therein) or condensate formation . The exper- 
imental results in Ref. W indicate, however, that as soon 



as the condensate becomes unstable the GPE alone is 
not suitable for the description of the dynamics of the 
collapse. In particular the observed substantial loss of 
condensate atoms contradicts the number conservation 
inherent in the GPE. This failure can be expected be- 
cause the collision term of the GPE does not account for 
collisional losses of condensate atoms due to non forward 
scattering |^ , which should become relevant for the rapid 
dynamics occurring during the collapse. 

More interestingly the gas was seen to be unstable even 
under conditions, where the GPE predicts the existence 
of a stable ground state. To clarify this aspect of the 
observations we have studied the stability of a Bose con- 
densed gas with a negative scattering length using the 
GPE, taking into account the role of the initial state, 
prepared in the experiment in Ref. We show that 
this initial ideal gas state leads to collapse condition that 
is different from that for the existence of a ground state 
of the condensate. 

Throughout this article the s wave scattering length of 
the interatomic force, oq, is assumed to be negative. For a 
Bose condensed gas in a spherical trap, with the harmonic 
potential Vtrap — ™-^ho''^/2i the natural units of time 
and length are the inverse trap frequency lo^^ and the 
harmonic oscillator length aho = \/^Jrrujj\^, respectively. 
Expressed through these units the GPE for an interacting 
Bose gas, consisting of N atoms, assumes the well known 
form |q| 
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Here the condensate wave function ■0 is normalized to 
unity, k is given by iV|ao|/aho and the corresponding unit 
of the energy is then Nh^ /ma\^. Equation (|l|) shows that 
k is the only free adjustable parameter of the GPE for a 
spherical trap. 

On first sight, the form of Eq. (|l|) might suggest that 
the critical conditions for collapse are governed solely by 
the parameter k. With respect to the existence of a stable 
ground state of a condensate this is true. In fact, the 
theoretical studies in Ref. PJlC|] confirm the condition 
for stability 



k <kcr^ 0.575. 



(2) 



When the initial condensate wave function deviates from 
the ground state, however, a collapse can occur even if 
the actual parameter k is well below the above standard 
value. 

As proposed by Pitaevskii j|] the conditions for stabil- 
ity or collapse can be visualized by means of two funda- 
mental estimates provided by the GPE. The first one is 
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the uncertainty relation between the mean square of the 
radius of a condensate, (r^) = J drr'^\ip\'^, and its kinetic 
energy, T = J dr\ xj ^\'^ /2, which reads 



(r^)T > 9/8, 



(3) 



where the choice of units corresponds to Eq. It is 
worth noting that once {r'^) vanishes the condensate wave 
function has coUapsed. The kinetic energy then becomes 
infinite. The second estimate consists in an inequahty 
between the mean field energy, E„-^{ = —Ank J dr\ip\'^/2, 
and the kinetic energy, 



(4) 



where /3 = 0.0575 is a universal constant. Equations (j^) 
and (^) are preserved for all times by the GPE. 




the mean square of the radius of a condensate and the 
trap energy are related through 



(6) 



The inequality given by Eq. then provides an estimate 
for the trap energy Eq. (^) through the kinetic energy. 
In addition, according to Eq. (Q), the mean field contri- 
bution to the total energy Eq. (|) can be estimated in 
terms of T. This in turn implies the general relation 



E > TH T" 
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which is valid for all times. 

In Fig. the lower bound Eq. (|^) is illustrated for the 
parameters k = 0.5 (Fig. |a), k = 0.53 (Fig. |l|b) and 
k = 0.575 (Fig. |l|c) and an initial ideal gas ground state. 
For this Gaussian initial condensate wave function the 
total energy and the initial kinetic energy are given by 
E ^ Eir, = 3/2 - fc/x/2^ and Tin = 3/4, respectively. 
In Fig. ^ the function /(T, k) exhibits a local minimum 
which encloses the point (Tin, E). As the total energy is a 
constant of motion and T{t) varies continuously in time 
the kinetic energy is caught in the valley for all times. 
As pointed out by Pitaevskii |^ , the inequality given by 
Eq. (^ then predicts the mean square of the radius of 
the condensate to be positive, i.e. (r^) > 9/8T > 0, and 
a total collapse, i.e. (r^) = 0, can never occur. As can 
be seen in Fig. |^c, very close to the standard critical 
interaction strength Eq. ^ the function /(T, k) becomes 
monotonous. Hence, a ground state of the GPE ceases to 
exist as soon as the local minimum of f{T, k) vanishes. 

In Fig. |l|b, however, the total energy is above the lo- 
cal maximum of /(T, k) and at least the lower bound 
in Eq. (0) allows the kinetic energy to access arbitrarily 
high values. In this case a total collapse of the conden- 
sate wave function cannot be excluded by the inequality 
Eq. (PI). 



FIG. 1. Lower estimate for the total energy i5 as a func- 
tion of the kinetic energy T for three different parameters 
of the interatomic force: (a) k = 0.5, (b) k = 0.53 and (c) 
k = 0.575. The horizontal lines indicate the total energy 
for an initial ideal gas ground state condensate wave func- 
tion. The vertical line indicates the initial kinetic energy of 
the Bose condensed gas. The energies are given in units of 
Nh^ /ma\^. The accessible points {T,E) are above the curve 
,f(T,fc). 

As an additional constant of motion, the total energy, 

E=.T + V + E^{, (5) 

contains the kinetic and the mean field part as well as 
the trap energy V — J drVtraplV'P- In the spherical case 
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FIG. 2. Numerical solutions of the time dependent Gross 
Pitaevskii equation for a spherical trap with an initial 
ideal gas ground state condensate wave function and dif- 
ferent parameters of the interatomic force: (a) k = 0.52 
and (b) k = 0.53. The graphs show the radial density 
p{r,t) — 47rr'^|-i/;(r, t)p as a function of the radius r and the 
time t. The time and the radius are given in units of the in- 
verse trap frequency lo^^ and the harmonic oscillator length 
Sho = \/ h/mu}\-ia, respectively. The radial density is normal- 
ized to unity. 

To clarify this situation, Fig. || shows the results of nu- 
merical simulations of the GPE in Eq. (Q) for fc = 0.52 
(Fig. ||a) and k = 0.53 (Fig. ^3), with an initial ideal 
gas ground state condensate wave function. In Fig. 0a, 



the simulation predicts the condensate to be stable with 
respect to collapse within the considered period of time. 
For the present parameter k = 0.52 the total energy is 
already slightly above the local maximum of the func- 
tion /(T, k) and a rigorous assertion concerning stability 
or collapse cannot be obtained directly from Pitaevskii's 
estimate. 
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FIG. 3. Comparison of the spatial condensate density 
p{r,t) = \ij){r,t)'^ at time t — 2.52 1^^^, close to the col- 
lapse, (solid line) with the initial Gaussian density (dashed 
line) for k — 0.53. The unit of the radius r is chosen as in 
Fig. ^ The normalized densities ( f drp{r,t) = 1) are shown 
on a logarithmic scale. 

In Fig. Ijb, however, the condensate wave function ex- 
hibits a collapse although the actual parameter of the 
interatomic force of fc = 0.53 is still sufSciently weak for 
a stable ground state of the GPE to exist. To visualize 
the contraction of the condensate wave function, in Fig. ^ 
a typical density |?/;(r, close to the time of the col- 
lapse, is compared with the initial Gaussian density. The 
corresponding wave function is obtained from the simula- 
tion in Fig. yb. As can be seen on the logarithmic scale, 
the predicted condensate density, close to the collapse, 
deviates from a Gaussian and exhibits a sharp peak at 
its center. 

The GPE thus predicts the preparation of an initial 
ideal gas ground state condensate wave function to lead 
to a pronounced reduction of the critical parameter of 
the interatomic force to 

fccr < 0.53 (8) 

in comparison with the standard critical value for the 
existence of a stable ground state. This critical value for 
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an initial ideal gas ground state has been reported also 
in Ref. [p^ as an intermediate result in the numerical 
study of stable ground states of the GPE. A numerical 
simulation alone, however, cannot predict stability of a 
condensate on arbitrarily long time scales. The required 
lower estimate of fccr is thus provided by the analysis 
illustrated in Fig. |l] and yields 

kcr > 0.5. (9) 

The sharply defined critical parameter of 

fccr = 0.461 ±0.012 ± 0.054, (10) 

reported in Ref. is indeed significantly smaller than 
the predicted standard value for the existence of a stable 
ground state in Eq. This corroborates the predictions 
of the GPE in view of the destabilizing role of the par- 
ticular initial state that was prepared in the experiment. 
Equations (H) and (^), however, are not strictly valid 
in the experimental case because the magnetic field was 
ramped linearly to its final value rather than abruptly. 
Interestingly, though, our result obtained from the GPE 
agrees with the observations within the experimental up- 
per limit for fccr. 

The analysis of stability or collapse of a dilute Bose 
condensed gas with a negative scattering length, per- 
formed in this work, can be applied to arbitrary ini- 
tial condensate wave functions and trap geometries. The 
optimal estimate from Eq. (0) is obtained for a spheri- 
cal trap. Related experiments with more general initial 
states should be feasible with the method described in 
Refs. and could provide a further understanding of 
the range of validity of the GPE in non equilibrium sit- 
uations. 
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